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Abstract 

We present the first measurement of tau-decays to hadronic final states with a </>-meson. This is 
based on 401.4 fb _1 of data accumulated at the Belle experiment. The branching ratio obtained is 
B(t± -» QK^v) = (4.06 ± 0.25 ± 0.26) x 10~ 5 . 

PACS numbers: 
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INTRODUCTION 



While hadronic decays with a <fi meson in the final state are valuable to investigate 
QCD at a low mass scale, they have never been observed due to their small branching 
fractions. The decay r ± — > ipK^v is Cabibbo suppressed and further restricted by its small 
phase space, while the decay r ± — > 0tt ± z^ is suppressed by the OZI rule although it is 
Cabibbo allowed. Taking into accounts differences in Cabibbo mixing effects and phase 
space relative to the Cabibbo allowed transition — > K°*K ± v, the branching fraction for 
t 1 * 1 — > $K ± v is estimated to be — > §K ± v) = 2 x 10" 5 On the other hand, the 
vector dominance model predicts Bij^ — > ^jtt^^) = (1.20 ± 0.48) x 10 -5 0. 

CLEO searched for these decays using 3.1 fb" 1 of data taken on T(4S) resonance. They 
set upper limits of B{t± -> 4>K±v) < (5.4 - 6.7) x 10" 5 and -> 0vr ± z/) < (1.2 - 

2.0) x 10 -4 at the 90% confidence level [l|. Here we report the first measurement of r 1 * 1 — > 
(j)K v and r — > (/>K 7r°u decays. These results are based on a data sample of 401.4 
fb _1 corresponding to 3.6 x 10 8 t + t~ pairs collected near the T(4S) resonance with the 
Belle detector at the KEKB asymmetric-energy e + e~ (3.5 on 8 GeV) collider 0|. The Belle 
detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector 
(SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov 
counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located inside a 
super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 0]. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
sample of 157.7 fb -1 , while a 1.5 cm beampipe, a 4-layer silicon detector and a small-cell 
inner drift chamber were used to record the remaining 243.7 1b- 1 |. 

While the decay r ± — > <pTT ± h / is also interesting, it is treated here as a background process 
including the kinematically allowed but phase-space suppressed decays of r ± — ► 0(n7r)^ 
(2 < n < 5). 



EVENT SELECTION 

We look for r ± — > ipK^v candidates in e + e _ — > t + t~ reaction in the following modes. 

7- s fgnai -»• + K± + (missing) 
^ K+K- 
T u. g -» We) T + n(< 1)7 + (missing) 

The detection of ^-mesons relies on the <ft —> K + K~ decay (B = (49.1 ±0.6)%); the final 
evaluation of the signal yield is carried out using the K + K~ invariant mass distribution. 

Selection criteria described below are determined, based on examinations of Monte-Carlo 
(MC) events. The background samples consist of rr (1570.0 fb -1 , which does not include any 
decay-mode with a meson) and other backgrounds such as qq continuum, £>°£?°, B + B~ 
and two-photon processes. For signal, we generate samples with 1 x 10 6 — > (pK^u, 
<j)K ± n i>, 07r ± z/ and (pTc^ir !/. 

Transverse momentum (p t ) for a charged track is required to be larger than 0.06 GeV/c 
in the barrel region (—0.6235 < cos# < 0.8332, where 6 is the polar angle opposite to 
the incident e + beam direction in the laboratory frame) and 0.1 GeV/c in the endcap 



4 



region (-0.8660 < cos# < -0.6235, and 0.8332 < cos9 < 0.9563). The energies of photon 
candidates are required to be larger than 0.1 GeV in both regions. 

To select r-pair samples, we require four charged tracks with zero net charge and a 
total energy of charged tracks and photons in the center-of-mass (CM) frame less than 11 
GeV. Furthermore, the missing momentum in the laboratory frame is required to be greater 
than 0.1 GeV/c and its direction to be within the detector acceptance, where the missing 
momentum is defined as the deficit of sum of the observed momentum vectors from that of 
the initial e + e~ system. We also require cos ^rust-miss < — 0.6 to reduce qq backgrounds, 
where ^thrust-miss * s an °P en i n g angle between the thrust axis and the missing momentum in 
CM frame. The event is subdivided into 3-prong and 1-prong hemispheres according to the 
thrust axis in the CM frame. These are referred to as the signal and tag side, respectively. 

In order to remove the dominant qq background, we require that the lepton probability 
Pn/e be greater than 0.1 and that the invariant mass of the particles on the tag side be less 
than 1.8 GeV/c 2 (~ m T ). Similarly, we require that both kaon daughters of the candidate 
have kaon probabilities Pk > 0.8. The effective mass of the signal side must be less than 1.8 
GeV/c 2 . Here Pi is the likelihood that a charged particle is of type t {t = /x or e or K or n), 
defined as Pe = Li/{L(, + L x ), where Lg and L x are the likelihoods of the particle for I and 
other species hypotheses, respectively, determined from responses of the relevant detectors. 
We allow at most one photon on the tag side to take account of initial state radiation, while 
requiring no extra photons on the signal side. 

Candidates mesons are reconstructed using oppositely charged kaons within the barrel 
and forward endcap region. To suppress combinatorial backgrounds from other r decays 
and qq processes, we require that the momentum be greater than 1.5 GeV/c in the CM 
frame. After these requirements the remaining contributions from B°B°, B + B~, Bhabha, 
\x pair and two photon backgrounds are negligible. 

To separate ipK ± v from 07r ± ^, the remaining charged track is required to satisfy the 
same kaon identification criteria as the daughters. The rr and qq contributions are 
reduced by requiring the opening angle (0™) between and in the CM frame to satisfy 
cos6™ > 0.92 and the momentum of the <j)K system in the CM frame must be greater 
than 3.5 GeV/c. For 07r ± z/, we require the charged track to be identified as a pion, P n > 0.8, 
and the opening angle between and 7r =1= in the CM frame to satisfy cos 8™ < 0.98. 

FigCJa) shows the K + K~ invariant mass distribution after all — > ipK^v selec- 
tion requirements. While there are two possible K + K~ combinatorial contributions from 
K ± K T K ± on the signal side in forming a meson; all combinations are included. Non- 
resonant backgrounds are mostly attributed to —>■ K + K^-k^v events with B = (1.55 ± 
0.07) x 10 -3 . Very small contributions are expected from qq processes. 

SIGNAL AND BACKGROUND EVALUATION 

The detection efficiencies e for — > (pK^u, 0ii' ± 7r t V and are evaluated, as listed in 
Table UJ from MC using KKMC |6j , where the V — A interaction is assigned on the vertices 
and the final hadrons decay according to non-resonant phase space. The signal — > <j)K ± i> 
detection efficiency is e^Ku = (1.82±0.01)%, including the branching fraction of — > K + K~ . 

Signal yields are evaluated by a fit to a p-wave Breit-Wigner (BW) distribution, convo- 
luted with a Gaussian function (of width a) to account for the detector resolution. The 
width is fixed to be = 4.26 MeV (PDG value 0), but a is allowed to float. From fits to 
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M(K,K) GeV/c 2 

FIG. 1: K + K~ invariant mass distributions for (a) — > ipK^v and (b) — * (pir^u. Points with 
error bars indicate the data. The hatched histograms show the expectations from t + t~ and qq 
background MCs. The open histogram is the signal MC with B(r ± — * (f)K^ 1 v)= 4 x 10~ 5 in (a) 
and B(r ± -> <\>tx ± v)= 6 x 1(T 5 in (b). The curves show the best fit results. See the text for details. 

TABLE I: Detection efficiencies e and cross- feed rates (%). 



Candidates 


Decay modes 
4>Kv (\mv §Kiftv 


T — > <f>Kv 
t — ► (frnv 


1.823±0.009 0.049± 0.002 0.327±0.006 
0.110±0.002 1.660± 0.014 0.009±0.001 



signal MC, the (f> mass detector resolution is found to be, a = 1.07 ±0.03 MeV/c 2 , which is 
consistent with the fit result to the data of 1.2 ± 0.3 MeV/c 2 . 

The K + K~ invariant mass distributions for — > (fiK^v and 07r ± z/ candidates are fitted 
with a p-wave BW distribution plus a linear and a second order polynomial background 
function, respectively, as seen in Figs. [Ha) and (b). The obtained signal yields are N T ^<pKu = 
573 ± 32 and N T ^ U = 753 ± 84. 

MC studies show that only r ± — > 07r ± z/, r ± — > (pK^TT !/ and qq samples yield significant 
contributions peaking at the mass. The contribution of other backgrounds is less than 
0.01% and can be neglected. The number of r ± — > (pTc^u events in the data is already 
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evaluated. Other contributions are estimated below. 

In order to evaluate the branching fraction and background contribution from — > 
<j)K ± TT°i>, we select 7r°(— > 77) candidates and combine them with ipK^v combinations that 
satisfy the requirements listed above. The signal yield is estimated by fitting its K + K~ 
invariant mass distribution with a p-wave BW distribution plus a linear background function, 
as shown in Fig|21 The resulting yield is 8.2 ±3.80i^7r°z/ events in data. With the detection 
efficiency of e^, Kn o v = (0.395 ± 0.007)% evaluated by MC and the produced number of rr's, 
N TT = 401.4 (ftr 1 ) x0.892 (nb) = 3.58 x 10 s , we obtain a branching fraction of 

Bfr* -> fiK^u) = (2.9 ± 1.3 ± 0.2) x 10~ 6 (1) 

with a systematic uncertainty of 6.9%, which is described later. Using this we estimate 
the contamination of r ± — > <pK ± 7r°v events in the r ± — > cf)K ± v signal as N^± 1 ^^ )K±7t0u = 
(6.8 ± 3.1) events, given the cross-feed rate of — > (pK^TT !/ to the r — > ipK^v signal to 
be (0.327 ± 0.006)%, as listed in Table □ 




1.03 1.04 1.05 
M(K,K) GeV/c 2 



FIG. 2: K + K~ invariant mass distributions for — > (pK^Tr^v. Points with error bars indicate 
the data. Histograms show the MC expectations of r-pairs (hatched) and signal (white) with a 
branching ratio of 3 x 10~ 6 . 

From the MC (731.1 fb _1 ) study, we find the contamination of qq is N q q = 12.1 ±4.5. To 
take into account the uncertainty in production in qq MC, we compare MC results with 
enriched qq data by demanding that the effective mass of the tag side be larger than 1.8 
GeV/c 2 . With this selection, the background is qq dominated and the other backgrounds 
are negligible. The yield in data is 262 ± 21 events and 117.4 ± 9.9 events for qq. We then 
scale the above estimate by the factor, / = 2.23 ± 0.47, and obtain the qq contamination, 
iV™ nt = 14.8 ± 3.5 events. 



RESULTS 

The peaking backgrounds described above, r — > (f)K ± Ti v and qq, are subtracted from 
the signal yield, N T ±^ K ± U = (573 ± 32) - (6.8 ± 3.1) - (14.8 ± 3.5) = 551.4 ± 32.3 events. 
To take into account cross-feed between — > ipK^v and — > t^-n^v due to particle 
misidentification [K <->■ ir), we solve the following simultaneous equations, 

N^Kv = 2N TT {e^Kv x B^k 

+ C x B<tmv ) , (2) 
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N^u = 2N TT 

{^4>Kv X B<t>Kv + £<tmv x -B</>7ri/) j (3) 

where B^ku and B^ v is the branching fraction for — >• <pK ± i> and — >• (p^u, re- 
spectively, e's are the detection efficiencies listed in Table |U is the efficiency for 

reconstructing r ± — > (pK ± v as r ± — > (pn^v while e^," is the efficiency for reconstructing 
— > 07r ± zy as r 1 * 1 — > cf)K ± v. The resulting branching fraction for r ± — > <pK ± i > is 

5^ ±l/ = (4.06 ± 0.25) x 10" 5 , (4) 

where the uncertainty is calculated with only statistical ones of N^kv an d The un- 

certainty in the detection efficiencies, e's, will be taken into account in the systematic error. 
B^ nu is obtained by the same way as 

B^ v = (6.07 ±0.71) x HT 5 , (5) 

but this is not the final branching fraction for the decay since the unknown contamination 
of t — > (p(mr)v {n < 5) decays still must be subtracted. 

Systematic uncertainties are estimated as follows: Evaluation uncertainties of the inte- 
grated luminosity, rr cross-section and trigger efficiency are 1.4%, 1.3% and 1.1%, respec- 
tively. Track finding efficiency has an uncertainty of 4.0%. Uncertainties in lepton and 
kaon identification efficiencies and fake rate are evaluated, respectively, to be 3.2% and 
3.1% by averaging the estimated uncertainties depending on momentum and polar angle of 
each charged track. To evaluate the systematic uncertainty of fixing in the BW fit, we 
calculate the change in the signal yield when is varied by ±0.05 MeV (the uncertainty 
quoted by PDG) The uncertainty is 0.2%. The branching ratio for (p — > K + K~ gives 
an uncertainty of 1.2% following the PDG 0. The backgrounds from N ( f >K ± u and N^ w ± u 
have uncertainties of 0.3% and 0.4%, respectively. The signal detection efficiency e^Kv has 
an uncertainty of 0.5%. A total systematic uncertainty of 6.5% is obtained by adding all 
uncertainties in quadrature. The resulting branching fraction is then 

Bir^ cpK^v) = (4.06 ± 0.25 ± 0.26) x 10~ 5 . (6) 

The systematic uncertainties for — > cpK ± T{ v are similar to those for r 1 * 1 — ► cpK ± v. The 
main differences are on the trigger efficiency of 0.3%, the detection efficiency e^Knu of 1.8% 
and the 7r° detection efficiency of 1.7%. Those provide a total systematic uncertainty of 
6.9%, and the branching fraction of 

Bfr* -> (pK±K°v) = (2.9 ± 1.3 ± 0.2) x 10~ 6 . (7) 

Finally, we examine a possible resonance state that intermediates the final (pK hadronic 
system. CLEO assumed a resonance having a mass of 1650 MeV/c 2 with a width of 100 
MeV/c 2 in the evaluation of their detection efficiency for r — > (pK^u, however no signal was 
found. We generate a resonant MC with the KKMC simulation program. The weak current 
is generated with a V — A form while the (pK system is assumed to be produced from a 
2-body decay of a resonance. In Fig. Efa), the (pK ± mass distribution for data is compared 
to MC; the combinatorial background is subtracted using the K + K~ sideband. Fig. H^b) 
shows the 0's angular distribution in the <f>K ± rest frame (cos a), where the momentum 
direction of (pK ± in the laboratory frame (P(<f> + K)) is taken as the reference axis. It 
indicates an isotropic distribution in the (pK^ system. For both the invariant mass and 
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angular distributions of <pK ± system, the phase space MC reproduces the signal distribution 
well. On the other hand, the 1650 MeV/c 2 state, indicated by the dotted histogram in 
Fig. Efa), clearly cannot account for the entire signal. Assuming resonant production, the 
best agreement with the data is found for a mass and a width of ~1570 MeV/c 2 and ~150 
MeV/c 2 , respectively, as shown by the dot-dashed histogram. However, since the shape of 
the resonant MC is similar to the phase space distributed MC, it is inappropriate to look for 
the intermediate resonance state with V ~ O(100MeV) in this narrow mass range of ~250 
MeV/c 2 . In fact, the phase space distribution (the open histogram) agrees well with data. 




FIG. 3: (a) invariant mass and (b) angular distributions for cpK system. The non-</>-resonant 
backgrounds are subtracted using the sideband spectra. Points with error bars indicate the data. 
Open histogram shows the phase space distributed signal MC, and dotted and dot-dashed his- 
tograms indicate the signal MC mediated by a resonance with M = 1650 MeV/c 2 and T = 100 
MeV/c 2 and M = 1570 MeV/c 2 and V = 150 MeV/c 2 , respectively. In MC, the branching ratio of 
4 x 10 -5 is assumed, (b) </>'s angular distribution in the (frK^ rest frame, where the direction of 
P{(f) + K) in the laboratory frame is taken as the reference axis. 



CONCLUSION 

Using a 401 fb _1 data sample, we report the first observation of the rare r decay mode, 
— ► (f)K ± i> ) with a branching fraction of 

-> (pK+v) = (4.06 ± 0.25 ± 0.26) x 10~ 5 . (8) 
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